Influence of modal loss on the quantum state generation via cross-Kerr nonlinearity 
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In this work we investigate an influence of decoherence effects on quantum states generated as 
a result of the cross-Kerr nonhnear interaction between two modes. For Markovian losses (both 
photon loss and dephasing), a region of parameters when losses still do not lead to destruction 
... of non-classicality is identified. We emphasize the difference in impact of losses in the process of 

state generation as opposed to those occurring in propagation channel. We show moreover, that 
. correlated losses in modern realizations of schemes of large cross-Kerr nonlinearity might lead to 

' enhancement of non-classicality. 
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T— I ] Nowadays cross-Kerr nonlinearity is considered as a promising tool for quantum computation and non-classical state 
C^l . generation [2]. An entanglement arising between modes participating in the cross-Kerr interaction can be used for 
' establishing an interface between matter qubits and 'flying' light photonic 01, for generation of non-Gaussian states 
, and Schrodinger-cat states !3] , and for performing quantum gate operations [H] ■ An interest to cross-Kerr nonlinear 
rS I interactions is heated up by both developing ways to implement effectively even very weak non-linearities (which one 
) \ ', commonly expects to have in practice) [1, |a, 0| i ^.nd by discovery of methods to produce sufficiently large self-Kerr 
^ ■ and cross-Kerr nonlinearities (such as implementation of the electromagnetically induced transparency (EIT) [1, [l^ 
^ ' and photonic crystals [3]). 

Decoherence is a main practical obstacle to implementations of schemes using Kerr and cross-Kerr nonlinearities. 
A genuine example of quantum state degradation due to losses is a decoherence of a quantum superposition state, the 
effect drastically enhanced if speaking of a macroscopic superposition state of the Schrodinger cat type. Already more 
• than 20 years ago the self-Kerr nonlinearity was proposed as a tool for generating such a Schrodinger-cat state (more 
[ precisely, a superposition of two coherent states with the same amplitudes but opposite phases) pjj]. However, photon 
losses turn this superposition into statistical mixture of two coherent states with the rate proportional to the square 
^vq modulus of the amplitude of these states. Modifications of the scheme for the cross-Kerr nonlinearities or four-wave 
^Kf^ i mixing brought no advantage with respect to photon loss 4]. This unfortunate circumstance made one look for the 
ways to circumvent the problem of decoherence that inevitably accompanies Kerr nonlinearity. Recent suggestions in 
this direction are based on the conditional preparation of desired states (which brings into consideration an additional 
problem of the finite detection efficiency), and are aimed to exploit weak nonlinearities [JjSQ- Recently, even a way 
[ to produce cat-states 'on demand' was suggested using a source of single-photons p^ . 
IL* ' In our work we want to discuss an aspect of the decoherence which has been seldom discussed when considering 
an influence of losses on states generated via Kerr nonlinearity. Namely, we address losses arising in the process 
of generation and not due to propagation of the generated state via lossy channels. We concentrate our attention 
5— i on a feature that might be quite signiflcantly pronounced in modern schemes of generating large Kerr nonlinearity: 
the modal loss can be strongly correlated. Indeed, the modes occupy the same volume and interact with the same 
physical systems which form the reservoirs. Also, if the Kerr-nonlinearity scheme implies a sufficiently strong dispersive 
coupling of light modes to emitters, then coupling of these emitters to dissipative reservoirs might also appear to be 
quite strong. As a result, this would mean strongly correlated modal losses. For example, in photonic crystals 
high density of states in the vicinity of a modal frequencies and emitter's transition frequency can cause the strong 
emitter-field coupling; but it would also imply higher population loss of emitters due to coupling to radiative reservoirs. 
Dephasing losses of emitters would as well invoke a correlated modal dephasing. 

Coupling to correlated reservoirs can drastically change state dynamics in comparison with loss to uncorrelated 
reservoirs. For example, it was demonstrated that coupling to the common reservoir preserves entanglement of a 
two-mode state [14]. Moreover, coupling to the common reservoir is capable of creating an entanglement between 
states of initially uncntangled modes even in absence of any direct interaction between them p^ . Il6j 

In our work we demonstrate both how the correlated loss arises via Kerr nonlinear process, and how it affects 
the generated states. For this purpose we derive analytic solution generalizing a powerful and illustrative method 
of Chaturvedi and Srinivasan [2fl|. On a number of examples we show how the correlated loss enhances and creates 
intermodal correlations and even entanglement, and might lead to generation of entangled states quite different from 
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those generated in the same scheme without loss. Correlated loss can result in the significantly enhanced robustness 
of the generation scheme. 

The outline of the paper is as follows. In Sec. II and in the related Appendices we describe how cross-correlation 
terms emerge via correlations of Markovian reservoirs; we consider an example of the emitter-field interaction schemes 
producing correlated modal losses in the Section III. Then in Sec. IV we describe the method for obtaining exact 
solutions of the cross-Kerr nonlinear interaction between modes in presence of losses to uncorrelated reservoirs and 
give generalization of the method for some cases of correlated losses. In the Sec. V we analyze infiuence of losses 
in the nonclassical state generation process for the case of uncorrelated loss. Some examples of correlated losses are 
considered in Sec. VI. 



II. MASTER EQUATION FOR CORRELATED AND UNCORRELATED LOSS 

To illustrate clearly an infiuence of correlated and uncorrelated losses, we restrict ourselves to the Markovian loss 
accounting for photon losses and dephasing of interacting modes. We start from the general effective Hamiltonian 
H(t) describing both self- and cross-interaction (for the moment we refrain from detailing it) and interaction of modes 
with reservoirs responsible for losses Vioss(0- 

V{t) = H{t) + Viossit), (1) 
Moss(i) = a\ri{t) + r\{t)ai + alr^it) + Tl{t)a2 + a{aiDi{t) + 4a2i?2(t). (2) 

Here we use the interaction picture with respect to the free Hamiltonians of reservoirs and the modes participating 
in the interaction process. These modes are described by usual bosonic creation and annihilation operators satisfying 

[ai, a|] — [a2, aj] — 1, [ai, = [0-1,0,2] ~ 0. 

The operators Ti 2{t) and Di 2{t) describe reservoirs responsible, correspondingly, to the photon losses in modes Oi 
and a2, and to dephasing of these modes. They may include also stochastic variables describing different realizations 
of reservoirs. 

It should be emphasized that reservoir operators ri_2(i) and £>i,2(0 (together with the initial state of the reservoir) 
completely describe the reservoir properties with respect to the interaction with the modes. These operators are built 
on the basis of underlying microscopic model and account for all relevant physical parameters. For example, if the 
photon loss reservoir of the first mode is composed of electromagnetic field modes with frequencies Wj , described by 
the creation and annihilation operators bj , bj , then 

= J^dj^i cxp{-iwjt}, (3) 
j 

where each gj is the constant of interaction of the mode oi with the jth mode of the reservoir. Sets of frequencies 
Wj and interaction constants gj describe completely physical properties of the reservoir. In particular, if the reservoir 
is the set of electromagnetic modes of a non-absorbing dielectric structure, they are found from the eigensolutions of 
Maxwell's equations for this structure [2^ . 

We make a number of standard assumptions in considering decoherence: coupling with the reservoirs is weak, the 
initial state of the modes Uj are uncorrelated with initial state of reservoirs, and correlation times of reservoirs are 
small enough to enable an implementation of the Born-Markov approximation. Then, using, for example, a time- 
convolutionless projection operator technique 18], one can obtain the master equation, the Liouville equation of the 
following form: 

-^p(0 = ~i[H{t),p{t)] + i[(7i - 7i2)/:(ai) + (72 - li2)Cia2) + (di - di2)Cia\ai) + (da - di2)/:(4«2)]p(t) + 

^[712^(01 -I- 02) + di2C{a\ai + ala2)]p{t), (4) 

where the superoperator C{b) acts on the density matrix as 

Cib)p{t) = 2bp{tp - bhp{t) - p{t)b'^b. (5) 

Here we set h = \ for simplicity. For more details of the derivation see Appendix A. By construction, the master 
equation ^ provides for non-negative definite p{t) for arbitrary t>Q. 
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Note that we call reservoirs "correlated", if their integrated cross-correlation function is non-zero, for example, 

and, correspondingly, the coefficient gi2 defined in Eq. (|A2[) is non-zero. Of course, one can always transform Eq. (|4]) 
to the diagonal form. However, in this case Lindblad operators (i.e. operators like b in the diagonal form ([5])) will 
be the linear superpositions of the former Lindblad operators, and the transformed equation will be still describing 
a coupling between physical objects represented by these original Lindblad operators (modes ai^2 in our case). For 
example, the possibility that both modes are coupled to the same reservoir (either the photon loss reservoir or the 
dephasing one) corresponds to an equality in relations (|A3p in Appendix A. In this case the operators describing the 
reservoir are proportional to each other, say, Ti{t) = xr2{t) and Di{t) — yD2{t). Equation ^ then reduces to 

^p{t) = -i[H{t), p{t)] + ^[jiC{ai + xa2) + diC{a\ai + ya\a2)]p{t). 

The most important point here is that in this equation both modes behave like a single object with respect to 
relaxation. Such a 'decoherence' is able to induce enduring entanglement between modes oi and 02 (an example is 
shown in the Appendix B). 



(a) (b) 




FIG. 1: Examples of level structures for generating cross-Kerr nonlinearities. Coupling to the quantized modes is shown by 
thin arrows and coupling to classical driving fields is shown by thick arrows. 

Furthermore, interacting light modes might experience correlated losses due to fact that they both interact with 
the same atoms. One can illustrate the mechanism of such a loss by the following qualitative consideration. Let us 
consider a general Hamiltonian describing the light-atom interaction Hq plus terms describing interaction of atoms 
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with the dissipative reservoirs 



^^totai — HQ{t; fli, 02, Sk) + ^ SjTj, 



(6) 



where Sj and Tj are operators describing the atoms and dissipative reservoirs, respectively. The Hamiltonian of this 
form describes a general interaction between field modes and emitters of some type as the correlated loss can occur 
in different physical systems. Provided atom-reservoir interactions in ^ is sufficiently weak as not to perturb much 
interaction between light modes and atoms, the following approximation can be used for the atomic operators: 



where U{t) 



S,it) « U^{t)S,{0)U{t) « Fj{t;ai,a2,Skm, 

Texp i^—i J dTHo{T; ai, a2, Sk)^ , and T denotes the time-ordering operator. After averaging out the 

atomic variable, the terms that describe reservoir-mode coupling in the effective interaction Hamiltonian will take the 
form of ^ fli, a2)rj. Here fj{t;ai,a2) — (Fj{t; ai, a2, Sk{0))) s,m and (. . .)am denotes the averaging over atomic 

j 

states (in general, over emitter's states). Thus, one can see that coupling of light modes to the same atom (emitter) 
interacting with the dissipative reservoir under the condition of adiabatic elimination of emitter's variables leads 
directly to the mode-reservoir interaction terms in the resulting effective Hamiltonian. Note, that these terms remain 
linear in reservoir operators Tj. Hence one can derive a master equation averaging over the reservoir in a standard 
way. In the Appendix B we give examples of the derivation of effective Hamiltonians and the corresponding master 
equations discussing the simplest two-mode Jaynes-Cummings system (Fig. HJa)). The described scenario of how 
correlated loss emerge is quite general, and can take place for a wide range of schemes involving light-shift-induced 
photonic nonlinearities. 

A three-level A-system interacting with classical driving and quantum fields represent a more realistic example of 
the scheme with the correlated loss. Consider the large cross-Kerr nonlinearity generation suggested in Ref. [lo| 
and depicted in Fig. [ijb). There two modes are coupled to the transition between 1 and 3 levels of the A-system 
in presence of two classical driving fields on transitions f-2 and 1-3. In the rotating- wave approximation and in 
the interaction picture with respect to free Hamiltonian of the reservoir, the problem is described by the following 
interaction Hamiltonian: 

Hctf = [gia\ exp{i(Ai - S)t} + g2a\ exp{i(Ai + 6)t} + Vl exp{iA2i} + ri3(t)]cri3 -I- Vl exp{iA2t}cr23 -I- H. c, (7) 

where £7 is the Rabi frequency of the driving fields; aui — \k) k,l = \, 2, 3; and 51^2 are the interaction constants for 
the coupling of a light mode to an emitter (atom) . Here for simplicity we have taken into account only losses on the 
transition 1-3. The setup depicted in Fig. [TJb) can be realized in toroidal microcavities, where oi and 02 correspond 
to the clockwise and counter-clockwise propagating modes [131 . As usually, we assume the Markovian reservoir and 
the following conditions hold: 

(ri3(i))r = 0, (ri3(i)rt3(T)), ^i5{t- r). 

From the Hamiltonian ([7]), the dynamics is governed by 

-j^CT3i « \[gia\ exp{i(Ai - 6)i} + g2al exp{i(Ai + 5)t} + exp{iA2t}](cr33 - an) + if7 exp{iA2t}CT2i. (8) 

Using the approach described in Ref. [lo| . we assume that level 3 remains practically unpopulated, (/^/(Ai ± (5) <C 1 
and Q,/ /S.2 <C 1, as well as gfc, 51 <C |Ai — A2I. Thus the A-system is prepared in the superposition of the metastable 
levels 1 and 2 (namely, in the state (|1) — |2))/\/2). Then, neglecting small and rapidly oscillating terms, one obtains 
from Eqs.QIHl) the following master equation 
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Ai-6 



a\ai + 



:4a2 ) ,p(<) 



Ai + ,5 



+ 2 [("^i ~ 712)^^(01) + (72 - 712)^^(02)] pit) + 



712 



£(ai -I- a2)p{t), 



(9) 



where 



71,2 = 7 



9I2 



(AiT<5)2' 



712 = 7 



5152 



(A?-<S2)- 
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The master equation ([9]) describes the cross-Kerr and self-Kerr interactions of two hght modes plus their couphng 
to the correlated reservoirs. As follows from the model depicted in Fig. mb)), it is actually the same reservoir: one 
can see that 7172 = 712- Naturally, loss rates 71.2 are much less than the loss rate of the emitter. However, if one 
deals with input modes in a coherent state of rather large amplitude for generating large cross-Kerr nonlinearities, 
then even comparatively small losses can strongly influence mode dynamics. Below we consider examples of such an 
influence, e.g., an example of a single- mode Schodinger-cat state. 

The occurrence of a correlated modal loss due to presence of the emitter (atom in our case) has been noticed in 
Ref. They have also pointed out that the loss rate should be proportional to the square of the ratio of the 

mode-emitter interaction constant and the detuning. However, we would like to emphasize that this loss rate does 
not depend on the population of level 3, as it can be seen from Eq. ([H]). 

III. CROSS-KERR INTERACTION MODEL 

Now let us turn to the specific nonlinear interaction described by the same type of the Hamiltonian as in the examples 
above and in the Appendix B. We will consider the effective Hamiltonian in the master equation ^ H{t) = Hq in 
the following general form 

2 

^^0 = ^ Xkia\aka\ai. (10) 
fc,;=i 

It describes the cross-Kerr and self-Kerr interaction with nonlinear coefficients Xki of two modes (or mode superpo- 
sitions) ai and 02. To solve the problem described by the master equation ^ with the Hamiltonian (jlOp . we adopt 
a simple and illustrative 'thermofield' notation [l9|, [20|. Essentially, instead of a density matrix acting on some space 
7i, say, p = "Yl, Pki\k){l\, where \k) is the Fock state with k photons in 7i, we consider a state vector \p) = 'Yl Pki\k)\l) 

k,l k,l 

in an extended space Ti^Ti* , where |/) is the Fock state with I photons in H* . So when the mode operators a and 
in the master equation ^ act on the density matrix from the left, one introduces operators a and is such a 
manner that 

\k){l\a > h^\k)\l), \k){l\a^ — > h\k)\l) . 

Obviously, operators a and commute with a and a) . An action of the superoperator C{a) on the density matrix 
can be represented in the thermofield notation as 

C{a)p — > \i{a)\p) — {2aa — a^a — a^a)\p). 

Also, the commutator of any function of the operators a and a^, e.g., the Hamiltonian HQ{t; a, a^) is re-written as: 

[Ho{t;a,a^),p] Ho {t)\p) = {Ho{t;a,a^) ~ Ho{t;a,d^)) \p). 

With help of these notations Eq. ^ can be represented in the 'Hamiltonian' form as 

4-Mt)) = Htotallp(i)) ^ 



dt' 



i(-2iHo(t) + (7ai - 7i2)L(ai) + (72 - 7i2)L(a2) + (^i - di2)L{a\ai) + {d^ - di2)L(4a2))|p(i)) 
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i(7i2L(ai + 02) + di2L{a\ai + a\a2))\p{t)) , . (11) 



Its solution is then of the form 



|p(t)) =exp{HtotaiOlp(0))- (12) 

The advantage of using thermofield notation over more traditional algebraic manipulation with superoperators is 
that in many situations (and, particularly, ones of our interest) it enables to simplify, make more illustrative and less 
cumbersome finding the solution (jl2p and estimation of time-dependent matrix elements. In particular, it allows to 
represent in a simple form a factorization of the superoperator exp{Htotaii} into multipliers with easily estimated 
actions on the number states 12011. 
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To illustrate this, let us consider a simple problem of modal loss in a single mode described by the equation 

The key to solving this equation lies in the observation that the operators 

A+ = a^d\ A- = aa, A3 = {a^ a + d + l)/2 (13) 

generate the SU(1,1) algebra with the Casimir invariant ^0 = a^a — a^a- Using the disentangling theorem [2lj for 
this group, we arrive at the simple result 

=exp{|ii}exp[-7,M3]exp[(l-e-^"*)A_]|p(0)). (14) 

Generally, the dynamics described by the solution (fT4| leads to transforming initially pure states into mixtures. 
However, for a coherent initial state with the amplitude a. 



\pm = \aw)= E exp{-iap} 

— „ y m\n\ 



'm) In) 



equation ITil) gives the following result: 



\p{t)) = |aexp{-7,V2})|a*exp{-7,i/2}). 

Returning now to the effective Hamiltonian (fTO|) that describes cross-Kerr and self-Kerr interaction of two modes (or 
modal superpositions) ai and 02, we write it in the thermofield notation as follows: 

2 

iJo ^ Ho = E XkiA^\2A^^ - 1). (15) 

Operators in Eq. (fTS]) are A^^'' = a^a^ — d\dk, ^3*"^ = (fllife + aj,.afe + l)/2. 

For completely uncorrelated reservoirs of different modes (i.e. 712 — di2 = 0) the master equation pT|) with the 
Hamiltonian Hq given by Eq. can be solved exactly using the following factorization [20] : 



exp{Htotaii} = exp 



dl ^4(l)^2 (a(2).2,- 4(1) 4(2) , 7l + 72 



X exp[(iPi ~ 7i)4^^^] exp[(iP2 - 72)4^^*] exp[ri_AL^^] exp[r2_ylL^^] (16) 



where superoperators A^'^^ are defined similarly as in Eq. (|13[) and 

Pfc = 2^XM, Tfc. = — — , Pk^2 2_^XkiAy. 

1=1 1=1 

It is useful to note that in Eq. all multipliers apart from two last ones are diagonal in the number-state basis. 
Also, it is easy to see that [F/j^, A^-*] = 0, operators T^^ are diagonal in the number-state basis, and operators 
are simply products of annihilation operators. Thus, Eq. (|16p provides for simple analytic solutions both for coherent 
initial states of interacting modes. 

The solution pH]) for uncorrelated reservoirs can be straightforwardly generalized for some special cases of correlated 
reservoirs and the Hamiltonian (fTO|) . For example, let us consider the problem without dephasing, dk — di2 — 0, and 
introduce the rotated mode operators bk as 

ai — hi cos((/)) -t- 62 sin(0), 02 — 62 cos(0) — bi sin(0), (17) 

where tan(20) = 2712/(72 — 71). Then for the non-unitary part of the master equation ()lll) one has 

(71 - 7i2)'C(ai) + (72 - 7i2)>C(a2) + 712-^(01 -I- 02) — > liCibi) + jC{b2), 

where 71 = 71 cos^(0) + 72 sm'^{(f>) — 712 sin(20), 72 = 72 cos'^{(f>) + 71 sin^((/)) + 712 sin{2(f>). Obviously, if the transfor- 
mation (|17p leaves the form of the Hamiltonian (jlOp invariant, one can derive an exact solution in the way described 
in this Section. 
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IV. UNCORRELATED RESERVOIRS 



A. General solution for the uncorrelated reservoirs 



In this Section we consider specific effects of tlie uncorreiated losses in tfie process of tfie cross-Kerr nonlinear 
interaction. Notably, the losses in such a nonlinear process can lead to loss-mediated correlations between the modes, 
as seen from Eq. ([16]) . These intermodal correlations modify the effect of losses on the quantum state generated in 
the cross-Kerr interaction with respect to the result, which one would intuitively expect treating the generation and 
loss separately, i.e. subjecting to loss a state that has been produced without losses. 

Consider a particular problem of generating an entangled two-mode state from initially uncorrelated coherent states. 
We assume that the nonlinearity is given purely by the cross-Kerr interaction, i.e. we put Xki = (1 — Ski)x/'^ in the 
Hamiltonian ([T0|) . Producing entangled states this way is important in a number of schemes of quantum computation 
and communication using continuous variables [3,0) HI. We assume modes oi and a2 to be initially in coherent states 
with amplitudes ai and a2, respectively. As was pointed in the previous Section, for this choice of initial states, the 
solution given by Eq. (jl6p has a simple form: 



\p(t)) = exp 



(2)^21 



f. ,/ t t s lit, I _i + 
exp < ixt(a[aia!2a2 — a[aia2a2) 7r((^i(^i + diO-i) 7r((^2'^2 + a2a2) 



X exp 



71 e- 



,ixt(a|ci2-a2''2)-7lt _ \ 



^Xiala2 - 4«2) - 71 



ai 



exp 



^2 ^QiXt(a{ai-a[ai)-j2t _ 



ix{a\ai - a{di) - 72 



a2 



|ai)|at)|«2)|«;). (18) 



A distinctive feature of the solution (|18p is that the exponential terms contain only operators diagonal in the number- 
state basis. Thus, in this basis the solution (fT8|) turns into 



where 



\p{t)) = e-l"^l 



-|Q2 



fc,/,m.n— 



„,„ , , exp<^ --[di(fc-0' + rf2(m-n)2]t 
\/k\l\m\n\ I ^ 



X exp ^ ixt{km -ln)-:f{k + l)~:f{m + n) + flll{t) + f^^\t) )\k)\l)\m)\fi) , 



Jmn \'^) 



ix{m - n) - 71 



72 



-I) -12 



\a2\ 



(19) 



(20) 



B. Analysis: When is the purity of the state not broken by damping? 

Now let us analyze the solution p9p in more detail and consider for the moment the case of no dephasing (di = 
d2 = 0)- In this case the expression for the purity of the state given by the solution (|19p is quite similar in structure 
to this solution itself: 



Tr{(p(t))n 



21. ^ g-2|ai|^-2|a2| 



E 



|ai|2(fe+')|a2|2("+") 



X exp 



{-7it(fc + - I2t{m + n) + 2Rc/(„2) (<) + 2Re/(;)(i)}. 



(21) 



To describe effects of simultaneous damping (photon loss) and cross-Kerr nonlinearity on the generated state, 
consider first a 'short-time regime' where the following relations hold for all the photon numbers k, I (m, n) that have 
a non-negligible probability of occurring in the state |ai) (|q;2)): 



\ixt{m - n) - 7it| < 1, \ixt{k -I)- i2t\ < 1. 
Then expressions (PU)) can be expanded as 

[ixt{m -n)- 7it]J (1) _ „ ,2 [ix^ (^ - - 72^]'' 



/i2)(t)=7itKPE' 
3=0 



U + 1)! 



/i:^(t)-72t|a2pE- 



j=0 



ij + 1)! 



(22) 



(23) 
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The first three terms of these expansions already describe quite well typical effects produced by simultaneous damping 
and cross-Kerr nonlinearity. For example, for fmn{t) in the limit of small times one has 



f!nlit)^7it\a,\ 



-, lit ,xt , \ , ^ \ 2,2 2,2/ \2i 1 ,2i \ 

1 - — +1— (to - n) + -Yl^t - X i (to - n) - -7iX^ (m - n) 
2 2 5 D 



Obviously, first three terms in the round brackets of this expressions do not lead to breaking of the purity of the state. 
Retaining only them in expansion renders unity value for the right-hand side of Eq. (|2ip . Indeed, assuming 



1 



^ » Ibft' - X't'im -nn 1 - ^ » Ibit' - X't'ik - l)% 



one obtains the time-dependent density matrix formally coinciding with the result for no photon loss 

«exp{-|ai(i)|n ["i(^)]^iW]' |fc)|[)|a2(t)e'^'-*)|a*(t)eix/^) 



k,l=0 



/kW. 



Here time-dependent amplitudes do not depend on the numbers k,l: 

Jit i 



ai{t) — ai exp 



72X|a2| O , a2it) = a2exp 



l2t i 



(24) 



(25) 



(26) 



The state given by Eqs. remains negligibly affected by losses if 71, 2^ ^ 1- Thus, the considered scheme of non- 
classical state generation is quite robust with respect to photon loss (in drastic difference with propagation losses of 
already generated cat-state where the off-diagonal terms |fc)(/|, A: 7^ I will decay with the rates proportional to 7|ap). 
In addition, it is interesting to note, that one might be able to satisfy conditions 71, 2^ <C 1 in schemes involving 
dispersive atom- field interactions in QED where it is possible to restrict losses to the photon loss of cavity modes ■ 
Moreover, further in this work we consider ways to circumvent an influence of losses by making them correlated. 

Remarkably, a purity of the generated bimodal state can be preserved not only in the case of small losses, but also 
for large loss. Indeed, in the limits of large losses one can consider x(fc — I) and x(to — n) as small quantities and 
expand functions (|20p in the following manner: 



\a2\-'fi]\t) « (1 - e-^^*) - ix{k - [(i + 
l«irVi'^(t)«(l-e-^^')-ix(TO-n) 



-72* 



72 
1 

71 



+ x^(to — n)^ 



t 



1 



72 72 

t 1 
- — + — 

71 7r 



-72* 



-71* 



1 



7? 



(27) 



Note that this approximation holds for arbitrary interaction times. If the interaction time is sufficiently large to fulfill 
the conditions 



l_e-72* »^2(^_;)2 

l-e"'^i* >X^(TO-n)^ 



1 



1 ^ ^ e-T2* ^ 



72 72 



1 



7l 



— ^ ^ ^ P-Ti* ^ 



71 7r 



7? 



(28) 
(29) 



then the state given by Eqs. (|19p is practically pure. Under the conditions (|29p this state is of the form described by 
Eq. (|25p with the time-dependent amplitudes given by 



o!i{t) = ai exp 
"2(0 = "2 exp 



2 

72^ 

2 



ix\a2\^ [ie-^^*-72"'(l- 
-ix|ai|Mte-"^*-7r'(l 



-72* 



-7lt\ 



(30) 



Clearly, purity of the resulting state in the long time-limit is precisely a consequence of a strong photon loss. In this 
way strong photon loss paradoxically suppresses state mixing predicted by the general solution (jl9p . 
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C. Survival of non-classicality for large losses 

There is another interesting feature that distinguishes the losses occurring in the process of Kerr interaction from 
the losses that take place after the interaction. In particular, the non-Gaussian state generated by the cross-Kerr 
interaction in the scheme discussed in Q can retain its non-classical features even for the loss level, which would 
completely eliminate any such features in case of free propagation of the state. To be specific, a typical signature 
of non-classicality of a quantum state is the fact that its Wigner function is negative in some regions of the phase 
space . We will show that for a 50% photon loss occurring in the scheme ^3)] during the cross-Kerr interaction, the 
Wigner function of the output state retains its negativity while the same photon loss occurring after the interaction 
would make the Wigner function necessarily positive. 




FIG. 2; (Color online) Wigner function W{u) of the non-classical non-Gaussian state of Ref. 0| exhibits strong negativity. 
Wigner function is viewed along the imaginary axis of the u plane. 



W(u) 









-2 2 4 6 

Re(u) 



FIG. 3: (Color online) Wigner function W{u) for about 50% photon loss in mode a corresponding to ^at = 0.7. The negative 
region of the Wigner function is now negligible, but still present, and the amplitude of the state is clearly damped. 



In the scheme Q for generating non-classical states, two coherent states in modes ai,a2 {a,b in notations of 
interact via cross-Kerr effect in a non-linear medium and subsequently the x-quadrature of mode 02 is measured. 
The resulting state of mode ai exhibits Wigner function with negative regions (see Fig. [J) and a characteristic 
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crescent (or banana) shape. Suppose mode oi is subject to losses during the cross-Kerr interaction, i.e., 71 > in 
Eq. while we assume 72 = c?i = 0^2 = 0. The photon loss in mode Oi is given by the reduction of the coherent 
amplitude described by Eq. ([26|) . Consider the situation when the mean photon loss is 50%. This corresponds to 
e~^'^* = 1/2 ^ 7it = ln2 « 0.69. The plot of the corresponding Wigner function is shown in Fig. [3] which shows 
clearly that although the negative region of the Wigner function is surpressed, it is still present. 

Now compare this with the situation when losses are introduced to mode ai after the lossless cross-Kerr interaction 
has taken place. Such losses are equivalent to mixing mode ai with the vacuum state on a beam splitter (BS) and 
discarding one BS output. A 50% loss corresponds to a 50/50 BS. It is known [23] that for such a balanced beam 
splitter the Wigner function of one BS output can be expressed as a scaled Husimi Q-function of the input state: 

VFout(a) =2Qin(\/2a). (31) 

The Q-function of a state p is defined as Q{a) — {a\p\a) /tt and is clearly non-negative for all coherent state amplitudes 
a, i.e., in the whole phase space of the mode. Therefore negative regions of the Wigner function cannot survive losses 
larger than 50% if these occur during propagation of the generated state. Hence losses that take place in the process 
of state generation via the cross-Kerr interaction are less harmful to the non-classicality of the output state than the 
same level of loss after the interaction. 



D. Dephasing into independent reservoirs 

As can be seen from the solution ()19p . an influence of dephasing into independent reservoirs can be profoundly 
destructive. Dephasing leads to diminishing of the off-diagonal elements in the number-state basis with rates propor- 
tional to the difference of these numbers. So for the coherent state with the amplitude Offc, a condition dfe|afc|^^ ^ 1 
should be fulfilled for the interaction time t to consider the influence of dephasing negligible. In the recently discussed 
QED schemes for generating non-linearity (including EIT-like ones), dephasing is usually disregarded without being 
estimated (see, for example, [l3|)- However, since the rate of losses increases with increasing intensity of the coherent 
states used in the generation process, more caution is required with respect to the dephasing. For schemes involving 
large cross-Kerr nonlinearity in the solid-state structures (such as, for example, photonic crystals) emitter-mediated 
dephasing could be a major source of state decoherence and a reason for failure of schemes involving initial coherent 
states with large number of photons. However, remarkably, if the cross-Kerr nonlinear interaction scheme is designed 
in such a way that losses due to dephasing are correlated, then it may be possible to avoid their destructive effects. 
This is the subject of the next section. 



V. CORRELATED RESERVOIRS 



It is well established, that the states of quantum systems can be correlated and even entangled through interaction 
with the common reservoir [l^, [l^ . This phenomenon can occur even in absence of any direct interaction between 
systems. In in Appendix B3 we give an example of such a phenomenon for a scheme of generating the cross-Kerr 
nonlinearity via dispersive interaction of modes with emitters. There a beam-splitting action of the common reservoir 
is considered, and it is shown how such a reservoir can produce a stationary entangled state of two modes. 

Here we focus our attention on another important possibility: namely, on a way to neutralize a destructive influence 
of losses in the process of generation by rendering these losses correlated and exploit a correlating effect of the reservoir. 

Let us consider an example of the realistic scheme to produce the large cross-Kerr nonlinearity described in Sectionllll 
(see Eq. (I?]) and the text thereafter). We consider the case of 

Ai-5 ~ Ai+S ^ ^^~K^' 
so the master equation ^ now transforms as 

+ \ [(71 - 7i2)/:(ai) + (72 - li2)C{a2)] Pit) + ^Cia, + a2)p{t). (32) 

In absence of decoherence (i.e., 7fc = 712 = 0), the scheme described by Eq. © is able to generate entangled 
superpositions of Schrodinger-cat states from initial coherent states of modes ai and 02 (for similar schemes see, for 



dt 



p{t) w ix 



[aXai+ala^^ , p{t) 



11 




FIG. 4: (Color online) Examples of the Wigner function for the conditioned cat state of the rotated mode bi. Figure (a) 
corresponds to absence of loss. Figure (b) corresponds to the perfectly correlated loss, 71 = 72 = 712 ~ lOx; for figure 
(d) 71 = 72 = 712 = 3x. Figure (c) corresponds to the completely uncorrelated loss, 71 = 72 = 3x, 712 = 0. Figure (f) 
corresponds to the partially correlated loss, 71 = 72 = 3x, 712 ~ 2.95x; figure (e) corresponds to the completely uncorrelated 
loss 71 = 72 = 0.5x, 712 = 0. For all figures x* = tt/S. 



example, [Tl|, [30|). It is easy to see that for xi = 7r/2 the scheme produces an entangled superposition of coherent 
states from a pair of initially uncorrelated coherent states, which reads as 



exp < 1 



I (a\ai + 402) I |ai)|a2) = (i|ai)|a2) + | - ai)| - 02)) 



(33) 



In the presence of correlated reservoirs one can find a solution of Eq. (|32|) noticing that under the rotation (|17|) the 
Hamiltonian part of Eq. (|32p remains invariant. Thus, performing the rotation one obtains 



i[7i£(&i)+72/:(&2)]p(t) 



(34) 



with the exact solution given by Eq. p6|l . Also, one immediately sees that for the completely correlated reservoirs 
(i.e. 71 72 = 712), the mode 61 is not affected by the loss, as 71 = then. Naturally, a cat state can be generated in 
this mode. Thus, in the limit of large loss, 72t ^ 1, and for xt = 7i'/2 it follows from Eq. that 



\p{7r/2x))^m^); I*) 



V2 



(i|Q;i)|Q:2) + I - ai)| - 5.2)) , 



(35) 
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where ai = aiCOS^((/)) — a2 cos((^) sin((/)), a2 = a2sin^(0) — ai cos((/)) sin((/)). We also assume that 2x\ai cos{(j)) — 
a2 sm(0)p ^ 72- From Eq. ([55)) . we derive a conclusion that the only effect of completely correlated loss is a reduction 
of amplitudes of the coherent states forming the superposition p3|) . 

Thus, we have seen that by making losses completely correlated one can completely avoid decoherence caused by 
these losses. Of course, in practice one can hardly have completely correlated reservoirs due to presence of additional 
uncorrclated loss (such as modal losses due to coupling to additional reservoirs etc.). Nevertheless, designing the 
scheme such as to have predominantly correlated losses might greatly enhance its robustness in production of non- 
classical states. We illustrate this with the simple example of the conditioned cat-state generation from the solution 
of Eq. ([9]). If the rotated mode &2 impinges on the detector, in the case of no signal on the detector the rotated mode 
bi is (up to the normalization factor) in the state: 

|pi(i))^ ^e^plixt{e-f)-:^ik + l) + Mt)\\k)\l), (36) 




where ai = ai cos{(j)) — sin((/)), ai — ot2 cos((/)) -I- ol\ sin(0), and 



7i(e'x*(fc-0-7i«-i) ^,(eix*(fc-0-72^ _ 1) 

/w(<) = TT = |ai| -I- — — laal . 

ix(fc - i) - 71 ix(k - - 72 

In Fig. 0] one can see examples of the Wigner function of the state (|36l) . For no loss (Fig. [D(a)) and = '''/2 
the state (|36)) is a usual Schrodinger cat state with the pronounced oscillations near the origin. Large correlated loss 
(71,2 ^ x) changes the size of the cat and rotates it (Fig. [H^b)), but otherwise leaves it intact. Lower correlated loss 
distorts the cat (Fig.HKd)) due to influence of additional mixing between modes in the interaction process (as it follows 
from Eq. ([55)1 ). Nevertheless, the state is strongly non-classical. Uncorrelated loss with the same rate eliminates the 
non-classicality outright (Fig. HJc)). Correlated loss allows for the non-classicality to survive (Fig. IH^f)) even if the 
uncorrelated loss with the rate equal to difference between individual rates and the correlation rate (e.g., 7 = 71 —712) 
destroys the non-classicality completely (Fig. UJ^e)). 



VI. CONCLUSIONS 



In what presented here, we followed the quest to find ways to impair the decoherence processes in quantum state 
generation and manipulation. For the particular class of nonlinear interaction processes, we have found two striking 
examples of loss dynamics, for which the losses themselves counteract decoherence: loss-mediated correlations between 
the interacting modes and losses to the correlated reservoirs. The latter result in strongly correlated modal loss. 

These correlated losses influence the dynamics of the modes undergoing the cross-Kerr nonlinear interaction in a 
completely different way than losses into independent reservoirs, the aspect of quantum nonlinear dynamics to large 
extent unexplored so far. Thus, remarkably, designing the schemes for the generation of the Kerr nonlinearity in such 
a way that losses in this process are correlated, one can greatly diminish their destructive impact and even exploit 
them for entanglement generation (see also Appendix B). As to the origin of this effect, if both modes, for example, 
interact with the same emitter transition, emitter losses are likely to lead to the correlated loss of both of these modes. 
Note, that in this case a significant modal loss might occur even in the case when the emitter subject to losses stays 
in superposition of metastable levels with negligibly small probability to occupy higher, decaying levels. For more 
details and examples on the origin of the correlated losses and their entangling effect see Appendix B. 

Turning to the other aforementioned unexpected aspect of the quantum dynamics, we have demonstrated that 
losses in the nonlinear process of state generation affect the quantum state in quite a different way to the propagation 
losses of already generated nonclassical state. This is mainly due to the loss-mediated correlations between the modes 
participating in the cross-Kerr interaction. In addition, losses through coupling to the correlated reservoirs (correlated 
loss) further enhance the difference between the decoherence processes during and after the state generation. In 
particular, non-classicality seems to be more robust with respect to the generation loss than to the propagation loss. 
We discussed an example of generation loss exceeding 50% with negative values of the Wigner function preserved, 
whereas the propagation loss exceeding 50% renders the Wigner function completely positive. 
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APPENDIX A: DERIVATION OF THE MASTER EQUATION 

We start from the general effective Hamiltonian H{t) of Eq. ([l] ^ describing both self- and cross-interaction and 
interaction of modes with reservoirs responsible for losses (see Section HI)) : 

V{t) = H{t) + MossW, Moss (i) = a{ri{t) + T{{t)ai + 4r2(i) + rl{t)a2 + a\aiDi{t) + 4a2i^2(t). 

Setting h — 1 for sake of simplicity, and using a time-convolutionless projection operator technique, one obtains the 
following master equation p^ : 

^p(i) = -i[H{t),pit)]- 
t 

j dr {(14oss(0l^loss(r))rP(0 + p(t)(Moss(r)l1oss(i))r - (T4oss(T)p(i)14oss(0)r " CHoss (t)p(O ^loss (T))r} (Al) 


where p{t) denotes the density matrix of the system averaged over states of all reservoirs (and over possible stochastic 
variables, too), and (. . .)]■ denotes averaging over all reservoirs. Reservoir modes are assumed to be initially in the vac- 
uum states. Further, we consider reservoirs of different types to be independent; correlation functions (T j{t)D\,{T)).c^ 
j,k — 1,2 are taken to be zero. We consider the reservoirs of the same kind to be mutually correlated, i.e., for 
Markovian dephasing reservoirs we assume 

(A(i)A(T))r = id,<5(i-r), {Diit)D2iT)), = idi2 5(t-r). 

Here S(t — t) is the delta-function and the rates di,d2, are real and non- negative. The cross-correlation parameter 
di2 is taken to be real. For the photon loss reservoir, we assume that the term a[Ti{t) + T\{t)ai + ajr2(t) -I- T2{t)a2 
in the Hamiltonian ([T]) preserves the total number of photons, i.e. only non-zero correlation functions are 

(r.(t)rl(r)), = S{t - r), (ri(i)4(r)>, = ^712 S{t - r). (A2) 

Also here, for simplicity, the rates 71 , 72 are assumed to be real and non-negative and the cross-correlation parameter 
712 to be real. For self-loss rates, 7^, di, and 'cross'-loss rates, 712, du, the following relations hold 

dadb > dlf,, -ialb > lib- (A3) 
Under the assumptions made above, the master equation of Eq. ^ is obtained from Eq. (jAip . 



APPENDIX B: ORIGIN OF CORRELATED LOSS 



1. An example: dispersive two-mode Jaynes-Cummings model with damping 

We will illustrate the process of appearance of a correlated modal photon loss in a process of off-resonant interaction 
between a mode and an emitter with a simple example. Consider two modes of the same frequency interacting off- 
resonantly with just a single two-level system (TLS), see Fig. [T](a). In the rotating- wave approximation and in the 
interaction picture with respect to the free Hamiltonian of the reservoir, in the frame rotating with the TLS transition 
frequency ujq, one has the following Hamiltonian describing the problem: 

HjK = A(4ai + 4a2) + (^^(ffiai + 52a2) + (ffiaj + 524)^") + (f^^r(t) + T+{t)a-) . (Bl) 

Here A = ujq — w and w is the mode frequency; gi^2 are interaction constants for the corresponding modes; and tr^ 
are Pauli operators for the TLS, cr"*" = |2)(1|, a~ = |1)(2|, az = (t^(J~ —a~a^\ vectors |fc), k = 1,2 describe the lower 
and the upper TLS levels, correspondingly. The reservoir operator r{t) describes the TLS energy loss. We assume 
this reservoir to be Markovian and the following relations hold: 



(r(i)), = o, (r(t)rt(r)), = 7<5(t-r) 
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where brackets (. . .)r denote an averaging over the reservoir. Here we are assuming that losses are weak (i.e., 7 ^ 

lffl,2|). 

We adopt the usual conditions for an adiabatic elimination of the emitter, i.e. the TLS starts at the lower level, 
and the detuning A between the mode frequency and the TLS transition frequency is much larger than 31,2- Thus, 
the TLS upper level remains practically unpopulated. Changing to the interaction picture with respect to the part 
of the Hamiltonian (IBip corresponding to the absence of the TLS-field interaction, we get the following interaction 
Hamiltonian 

VjK{t) = G (CT+Cexp{iAi} + h.c) + (fT+r(t) + H. c.) (B2) 
with the bosonic annihilation operator for the collective mode 

C = I } {giai + 3202) 

and G — \/ g\ + (7I • A formal solution for a^(t) without losses up to the third-order terms can be approximated as 

t 

a+(i) « a+(0) + i / Ati{2F\tr)F(,t{) - (B3) 



G 

where X{t) ^ GC(t) exp{-iAi} and F\t) w cr+(0) + cr3(0)C^(t) — (1 - exp{-iAO), if one takes into account the fact 

that the modal dynamics is very slow on the scale of the TLS dynamics and crz(t) can be considered as practically 
constant. From Eq. (jB3p . neglecting small and rapidly oscillating terms, after averaging over the atomic variables 
one arrives to the following effective interaction Hamiltonian 

VjK{i) ~ 2^CHt)Cit) + i^cHt)Cit)CHt)Cit) + ^ [CHt)T{t){l - exp{~iAt}) + H. c] . (B4) 

Deriving the master equation in the standard manner, one obtains an equation describing the correlated photon losses 

^p(t) « -i[SwC^C + xC^CC^C, p{t)] + jC{C)p{t), (B5) 

where Sw = 2G^/A, x = 4GVA3 and 7 = 2-fG^/A^. 

Effectively, Eq. (|B5|) describe both coupling between modes and their interaction with the same reservoir. Both 
these interaction might lead to the entanglement between modes. As seen in Section V, even in the absence of direct 
intermodal coupling (i.e. for Sw = 0, x = 0) an interaction of modes with the reservoir entangles these modes. 

The analysis made above can be readily generalized to different schemes of cross-Kerr nonlinearity generation 
through interaction of two modes with the same ensemble of emitters [l,[l3|. In Subsection |TT] we devise the procedure 
for the scheme of the giant cross-Kerr nonlinearity generation suggested in Ref. p^ . 



2. An example: dispersive two-mode Jaynes-Cummings model with dephasing 

Here we illustrate an appearance of correlated modal dephasing with the example of the Jaynes-Cummings model 
(Fig. [T|) considered in the previous Subsection. We model an influence of the dephasing reservoir as a stochastic 
fluctuation of the TLS transition frequency. In the rotating-wave approximation the problem is described by the 
following Hamiltonian: 



HjK = w{a\ai + 0^02) + i[cc)o + C{t)]'^z + cr^(.9iai + .g2a2) + {gia\ + g2a\)<7 



(B6) 



where C,{t) is random process describing a small rapid stochastic modulation due to non-radiative interaction with 
surroundings. For simplicity we take C,{t) to be just a white noise satisfying the following relations 



(c(t))s = o, {amr))s = d5{t~T), 



where (. . .)s denotes classical averaging. We consider the case of the weak loss, d <C |ffi,2|- 
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As before, we assume that the conditions for adiabatic elimination of the emitter hold. In the interaction picture 
with respect to the part of the Hamiltonian (|B6p corresponding to the absence of the TLS-field interaction, we have 
the following interaction Hamiltonian 

VjK{t) = G {<7+Cf{t) + H. c.) , fit) = exp + i / "^^^^^^ | ' (^^) 

A formal solution for (7^{t) in this interaction picture can be approximated as 

t 

a+{t) « a+(0) +iGj dh[2F^{h)F{h) - IjCHh)/* (h), (B8) 



where 

F+{t) « a+(0) + iGaMCHt) t dT/*(T). (B9) 



After averaging over TLS states, the following effective interaction Hamiltonian can be obtained from Eq. (|B9[) : 

VjK{t) « 2G^C\t)C{t)p^^\t) + A^C\t)C{t)CHt)C{t), (BIO) 

where 



P 



(2)(t)=Re 



f{t) / dT/*(T) 







Averaging over dephasing noise and neglecting small terms, we get the following master equation using the standard 
technique implemented to derive Eq. 



where 



-^pit) « -i[(5wC^C + xC^CG^C, pit)] + dC{C'^C)p{t), (Bll) 
at 



d = 4G^{ I drp{Q)p{r)\ ^ d^. (B12) 



The calculations of the coefficients in (|B12p are carried out using the following property [1^, [2; 

(■ 



t T 

dr J da;(C(r)C(x)), ^ =exp{-dt} 





and the fact that the detuning A is assumed to be large, A ^ d. 

So, one can see that dephasing of the atom leads to appearance of the correlated modal dephasing practically in 
the same manner as atomic population losses lead to the correlated modal loss considered in the previous Subsection. 
Also, modal dephasing occurs notwithstanding the fact that the upper atomic level remains practically unpopulated. 



3. Beam-splitting by decoherence 

Finally, we demonstrate that the scheme described by the master equation jBS} can effectively produce entanglement 
between the modes. In fact, this scheme can act as a kind of lossy beam-splitter even in the absence of intermodal 
interaction in unitary part of Eq. (jB5p . Indeed, let us consider a completely uncorrelated single-photon initial state 
|5'(0)) = |l)i|0)2, and the initial density matrix |p(0)) = |^'(0)) |^(0)) . In the zero- and single-photon subspaces one 
can assume the following orthonormal basis: 

= ^(5l|l>l|0)2 +52|0)l|l)2), IV--) = ;!^(52|l)i|0)2 - 5l|0)l|l)2), \v) = |0)l|0)2. 
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FIG. 5: Examples of the negativity, N{p), dynamics given by the solution (|B13[) for the initially disentangled state of modes 
1 and 2 (namely, single photon in the mode 1 and vacuum in the mode 2). The time, T, is given in units of (72; solid, dotted 
and dashed lines correspond to 7 = 32, 0.25(;2, 2^2 and gi ~ g2- Dash-dotted line corresponds to 7 = 32 and gi = 2g2. Ifere 
Slu + X = for all graphs. 



One can easily see that the state j-i/'-) is not affected by the losses described by Eq. (jB5|) because C|^_) = 0. Also, 
the following relations are satisfied: 

Thus the system of equations for the density matrix elements can be obtained from Eq. (|B5[) : 

W = P— (0), = -27P++(i), ^P+-W = -[7 + i(^w' + x)]p+-W, 

p_,(0=P-.(0), -^n+,{t) = -[^ + ii5w + x)]p+vit). (B13) 
at 

The solution (|B13p describes an emergence of entanglement form the initially uncorrelated state of both modes ( 
single photon in the mode 1 and vacuum of the mode 2). Figure [5] depicts a measure of entanglement, a negativity 
as given in Ref. [28.] 

7V(p) = i(TrV^-l), 

where a is the density matrix p partially transposed with respect to the first mode. Non-zero value of the negativity 
means that the state is the entangled one. It can be seen that the decay rate into the common reservoir does not affect 
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the finally reached entanglement. This rate affect only time during which a stationary state is reached. It follows 
from this system of equations that the initial uncorrelated state |l)i|0)2 under the action of the correlated modal loss 
asymptotically turns into 



The state (|B14[) is entangled for an arbitrary 2-^0- However, the maximal degree of asymptotic entanglement is 
reached when gi = g2, and with increasing of difference between gi and g2 the asymptotic entanglement decreases 
(Figure El). 

Note that the state (jB14p is influenced neither by the cross-Kerr interaction of the modes nor by the linear excitation 
exchange, and one can set both Sw = x = 0. The same type of state is produced by a correlated dephasing described 
by Eq. (|B11|) . Effectively, the correlation of reservoirs allows for existence of decoherence-free subspaces to which the 
two- mode state eventually evolves [l5j | . Entangling through the common reservoir with appearance of the long- living 
state similar to the one described by Eq. (|B14p might occur for emitters and collective reservoir modes near the 
band-edge in photonic crystals [29^ . 
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